Spin-blockade qubit in a superconducting junction 



C. Padurariu^ and Yu. V. Nazarov^ 

^ Kavli Institute of NanoScience, Delft University of Technology - Lorentzweg 1, 2628 CJ, Delft, The Netherlands 



(N 

O 
(N 

O 

O 

oo 



PACS 74.45.+C - Andreev reflection (superconductivity). 

PACS 74.50.+r - Josephson effect, tunnelling phenomena (superconductivity). 
PACS 73. 63. -b - Electronic transport in nanostructures. 

Abstract ~ We interpret a recent pioneering experiment [Zgirski M. et al, Phys. Rev. Lett., 106 
(2011) 257003] on quasiparticle manipulation in a superconducting break junction in terms of spin 
blockade drawing analogy with spin qubits. We propose a novel qubit design that exploits the 
spin state of two trapped quasiparticles. We detail the coherent control of all four spin states by 
resonant quantum manipulation and compute the corresponding Rabi frequencies. The read-out 
technique is based on the spin-blockade that inhibits quasiparticle recombination in triplet states. 
We provide extensive microscopic estimations of the parameters of our model. 
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Introduction. — The spin degree of freedom pro- 
vides a natural representation of quantum information in- 
spiring some of the first qubit designs Nanodevices 
that realize read-out and manipulation of single electron 
spins are in the focus of modern research [2]. An alter- 
native is provided by Josephson-based superconducting 
qubits [sj where the qubit states emerge from the inter- 
play of Josephson effect and Coulomb blockade, not hav- 
ng anything to do with real spin. There are theoretical 
proposals [4]|5j attempting to combine the advantages of 
both qubit realizations. This can be achieved by a su- 
perconducting spin qubit, where quantum information is 
stored in the state of a spin trapped in a superconducting 
junction. 

Generally, in this situation one deals with the spin of 
a superconducting quasiparticle rather than with an elec- 
tron spin. The supercondcuting quasiparticles are exci- 
tations of the Bardeen-Cooper-Schrieffer (BCS) supercon- 
ducting ground state that are spin doublets despite the 
fact that they do not bear a definite charge, being super- 
positions of electron and hole excitations (6l. 

One needs to localize a quasiparticle to control its spin 
state. It is possible to trap a quasiparticle not only inside 
a superconducting island [t], but also in eventually any 
superconducting junction. The quasiparticles are kept in 
the Andreev bound states (ABS) that develop [s] in the 
presence of a superconducting phase difference dropping 
at the junction. The presence of an additional quasipar- 
ticle in a bound state is manifested by a change of the 
superconducting current. In addition, spin-orbit interac- 



tion in the junction makes the superconducting current 
sensitive to the spin state of the trapped quasiparticle [4] . 
This provides a way to detect the quasiparticle occupancy, 
as well as to manipulate the spin state. Using ABS is par- 
ticularly advantageous since in this case the junction may 
support larger superconducting current as compared to 
the systems involving Coulomb islands [9]. 

Controllable trapping of quasiparticles in ABS and their 
detection has proven an experimental challenge and has 
been realized only recently 10 . The experiment uses 



a mechanically controlled break of Al strip to produce 
atomic-size junction of adjustable transparency. The junc- 
tion normal state conductance is contributed by a few 
transport channels; their transparencies are accurately 
characterized from the fit of TV curves. In the super- 
conducting state, each channel gives rise to an ABS that 
can accommodate quasiparticles. 

Detection is based on the fact that the critical current of 
the junction is altered by the presence of trapped quasi- 
particles. The measurement of critical current must be 
achieved with a time resolution smaller than the quasipar- 
ticle escape time. The experiment uses rectangular pulses 
of bias current of typical duration ~ 1 /is and different 
amplitudes. Each pulse induces with a certain probabil- 
ity the switching of the junction from the zero-voltage to 
running state. The latter is manifested as a measurable 
voltage response. The probability is determined by re- 
peating the measurement a big number of times ~ 10** for 
each value of the pulse amplitude [10| . Upon increasing 
the pulse amplitude the switching probability increases in 
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steps. The step positions indicate the critical current of 
the junction with different numbers of trapped quasiparti- 
cles and are in good agreement with the theoretical predic- 
tions. In addition, the measurement provides information 
about the relative probability of different occupancies of 
the ABS. These probabilities are proportional to the in- 
crease in switching probability at the corresponding step. 

The following experimental detail has inspired the pro- 
posal presented in this Letter. In most cases either zero 
or one quasiparticles (at lowest level) have been detected. 
However, the state with two trapped quasiparticles has 
been reported as well. It could not be observed in a junc- 
tion with a single high-transparency channel, but has been 
observed for junctions supporting two open channels. Im- 
portantly, the measured critical current corresponded to 
the case when the quasiparticles are situated in different 
levels: the lowest one and the next-to lowest one. 

In this Letter, we interpret this detail as the spin- 
blockade of two quasiparticles. Spin-blockade in a dou- 



ble quantum dot connected to two leads (11-13 is utilized 



for the read-out and operation of the most practical spin 
qubits. The transport cycle through the dots is arranged 
in such a way that the state with two electrons in differ- 
ent dots can only transit to a spin-singlet state with two 
electrons in the same level in one of the dots. Spin con- 
servation forbids transitions between singlet and triplet. 
Therefore the electron transport is blocked if two electrons 
in two different dots happen to be in a triplet state. 

Similarly, two quasiparticles in two different Andreev 
levels : lowest and next-to-lowest — form the ground state 
of the system with full spin S ~ 1. Spin conservation in 
this case forbids a transition to S — no-quasiparticle 
state of the lower energy, the transition proceeds quickly 
if the two-quasiparticle state is a singlet. One can say 
the junction is blocked in a triplet state. This guarantees 
long life-time and therefore observability of these specific 
two-quasiparticle states. 

We show that this spin blockade in combination with 
spin manipulation techniques can be used to control and 
read out a simple qubit. We develop a minimum the- 
oretical model that includes all relevant spin-dependent 
phenomena in the junction and provide extensive micro- 
scopic estimations for the parameters of the model. We 
detail the resonant quantum manipulation of all four spin 
states formed by two quasiparticles and outline a "natu- 
ral" read-out in the system. Experimental realization of 
our proposals would for the first time unambiguously prove 
the spin properties of superconducting quasiparticles and 
open up new ways to combine spin and superconductivity 
in the context of quantum technologies. 

Theoretical modeL — The break junction is char- 
acterized by a set of transport channels labelled by p 
with transmission coefficients < Tp < 1. The An- 
dreev levels corresponding to the channels have energies 
Ep{ip) ^ A(l - rpSin2(v3/2))i/2, where A is the bulk su- 
perconducting gap (~ 1.9 X 10"'* eV in 10 ). Each level 




Fig. 1: Andreev bound states in a junction with two open 
channels, a) Sketch of the junction characterized by the chan- 
nel transmission coefficients {Ti,r2} and the superconducting 
phase difference f. b) Energy of Andreev states as a function 
of ip. The three plots show most long-lived states stable with 
respect to spin-conserving energy relaxation. They are: the 
state with no quasiparticles (5* = 0), the state with a single 
quasiparticle (S = 1/2), and the state with two quasiparticles 
in a triplet spin configuration (5* = 1). c) The energies of the 
four spin states of two quasiparticles plotted versus the spin 
magnetic field. Inset: Zoom in the low _B-field region to show 
the spin-orbit splitting of the triplet states. The energy of the 
spin-orbit splitting is small compared to the exchange splitting 
esc ^ J ■ 



may accommodate up to two quasiparticles of opposite 
spin. We introduce the level occupation rip = 0, 1, 2. 

The occupation of the Andreev levels determines the 
phase-dependent part of the junction energy, which is 
given by the sum over the levels. 



i?((p)=5](np-l)£;pM, 



(1) 



and the superconducting current flowing in the junction 
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We focus on the case when two lowest levels are occupied 
with a single quasiparticle each, ni = n2 = 1, rip = if 
p ^ 1,2. We call these states TQTL (two quasiparticles, 
two levels). In this case, the two lowest levels do not 
contribute to the current resulting in a suppression of its 
magnitude. Two spin-1/2 quasiparticles in a TQTL give 
rise to four spin states: a singlet \S) and three triplet 
states \Tj), j = -1,0, +1. 

Since the TQTL is an excited state , one needs to drive 



the junction out of equilibrium to realize it. In 10 this has 



been achieved by applying a current pulse of large ampli- 
tude just before each measurement of the critical current. 
The current in the pulse exceeds a critical one and re- 
sults in quasiparticle generation at the junction. After the 
pulse, some of the quasiparticles generated are trapped in 
the ABS. Two quasiparticles are trapped with some prob- 
ability P2- The quasiparticles can relax their energy by 
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going to the lower Andreev levels. If their state is a sin- 
glet, they may also annihilate: we assume that this process 
is fast at the scale of the measurement time. The anni- 
hilation takes place also in the case of a singlet TQTL: 
in this case, we characterize it with the rate F^. If the 
quasiparticles are in the triplet state, they relax to TQTL 
and stay there manifesting themselves in a lower critical 
current. It is reasonable to assume random initial spin: 
in this case, the probability to measure the lower current 
is P = (3/4)P2- This is provided that the measurement 
time is shorter than the spin-orbit relaxation rate of the 
triplet state to be estimated below. 

In addition to the decay of the singlet state, the evolu- 
tion of the TQTL is determined by three spin-dependent 
effects: the exchange interaction, the spin-orbit (SO) cou- 
pling and the interaction with an external magnetic field 
B. We summarize these effects in a Hamiltonian Hs- 



exchange 



so 



(3) 



The exchange interaction term is expressed as usual in 
terms of the two spin operators cti and (J2 at the levels 1,2, 
-^exchange = ' f?2), J bciug the exchange coupling. It 
is convenient to rewrite the term in the basis of singlet 
and triplet states. 



^^cxchangc = - J/2 \S) {S\ + J/2^ |T,) (T, 



(4) 



The effect of SO interaction on a spin state of a quasi- 
particle occupying an Andreev level has been discussed 
previously in [4|[5]. The interaction is characterized by 
a level-specific pseudo-vector of dimension energy that 
polarizes the qusiparticle spin. 

We denote the SO pseudo-vectors of the two Andreev 
levels involved by ei and ^2 ■ The vectors generally depend 
on the superconducting phase drop ip over the junction. 
The following symmetry holds: ei.2{,^) = — ei,2(— f/?)- It is 
a consequence of time-reversal symmetry preserved by the 
SO interaction. The SO Hamiltonian takes the form [s] 
H'so — X]n=i 2 ^n(v) ' ^n- To cxprcss it in singlet-triplet 
basis it is convenient to introduce the sum and difference 
of pseudo-vectors e^^d — (ei ± e2)/2, following 
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Hso =Y.3's \Tj) {T,\ + ^ '-^^ \To) (Till 4 
^5 \S) (Tol + E ^^^^ + 



(5) 



The spin-quantization axis z is not yet fixed here. We will 
make a particular choice of the quantization axis further 
in the text. 

We also take into account the spin interaction with an 
external magnetic field = X]n=i 2 ^ ' '^n- H^re, the 
magneton /i is smaller than its value in vacuum. The 
reason for this is that the Andreev bound states may 
spread into superconducting leads at distance exceeding 




Fig. 2: Sketch of the atomic-size break junction. The wave- 
function of quasiparticles trapped in localized Andreev states 
extends away from the atomic-size constriction of length a to 
a much longer distance of the order of the superconducting co- 
herence length ^. Assuming the dirty-limit for the Al film su- 
perconductor used in the experiment, the motion of the quasi- 
particle inside the volume taken by the state follows a chaotic 
scattering trajectory. 



the screening length of the magnetic field. In what follows 
it is convenient to set jiB — >■ B, introducing a field B of 
dimension energy. In the singlet-triplet basis Hb takes the 
form 

Hji=Y,3B'\T,){T,\+ (6) 
y-^ B"" ± iBy , , , , , 

The magnetic field dependence of the energies of the 
resulting states is given in Fig. [TJ;. 

The dynamical evolution of the four spin states is ad- 
equately described by a density matrix equation which 
takes into account the unitary evolution described by Hs, 
as well as the decay of the singlet component with the 
rate Td- The density matrix p is defined in the space 
of four spin states of interest, pxy = {x\p\y), x,y G 
{S,T_i,To,T+i}. Owing to the decay, Trp < 1. The 
equation set reads: 



— p = i[Hs,p] 



TdPss\S){S\- (7) 



J2i^PsTAS) {T,\ + h.c. 



There are four equations to solve for diagonal components 
, as well as six for off-diagonal ones. It is relatively easy to 
obtain the numerical solution for arbitrary values of the 
parameters. To reveal the physics of the spin system we 
present analytical solutions in certain parameter regimes 
relevant for experiments similar to 10 



Before doing this, we need to estimate the parameters 
of our model to reveal the relevant time scales. 

Estimations. — Let us start with the energy relax- 
ation mechanisms in the junction. These mechanisms may 
involve the transitions of quasiparticles to lower Andreev 
levels, as well as their annihilation. The relaxation pro- 
cess dumps energy into an environment and thus must be 
accompanied by an emission of a photon, or of a phonon. 
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Since the states are locaHzed inside the superconductor 
where electromagnetic field is screened, the phonon mech- 
anism dominates. 



For the junction setup in 10 the situation is depicted 



in Fig. [2j The wavefunction of a localized quasiparti- 
cle extends in the superconducting leads for a length of 
the order of the superconducting coherence length, ^. As- 
suming the dirty- limit for the Al film superconductor [14| , 
the coherence length can be estimated as ^ ~ {hD/Ay^"^ 
where D is the diffusion coefhcient of electrons in nor- 
mal metal. The atomically-thin constriction extends for 
a length a ~ 1..10 x 10^^" m, therefore a ^ For the 
film geometry, the volume taken by the wavefunction of 
the ABS can be estimated as V ~ ^^d, d being the film 
thickness. 

Since the volume V is large at atomic scale, the rate 
can be estimated from the rate of phonon-mediated en- 
ergy relaxation for normal-metal electron states forming 
a continuous spectrum. This rate at energy E is given 
by Tr{E) — Kci-phE^, Kci-_ph being the electron-phonon 
coupling constant (Kci-ph — 3 ^s~^K~^ |15 for Al.) For 
our estimation, we need the transition rate Fjj between 
two discrete states. To get this, we divide F^ by number 
of states N available for relaxation, N ~ i^EV, v being 
the electron density of states. Therefore, F^ ~ T^/N . To 
adjust the estimation to the ABS case, we estimate E by 
A and V by the volume taken by an ABS. 

If we express ^ in terms of diffusion coefficient, the es- 
timation of N appears strikingly simple, N ~ {R\jG q)~^ , 
Gq = ire^ /h, R\j being the normal-state film resistance 
per square. One can argue that the same estimation holds 
for an arbitrary geometry provided R\j is replaced with 
an effective resistance i?cff of a metal piece of size ~ ^ ad- 
jacent to the constriction. This is a manifestation of the 



famous Thouless relation 19 between level spacing and 



dwell time in a piece of disordered metal. 

With this, we estimate the ABS relaxation rate as 

Trf — {GQRcs)Tr{A) ~ Kol-phA'^ GgiJeff • 



(8) 



This estimation of the phonon relaxation rate shall be 
compared with the results of Ref. [Tt] where F^ ~ F^ (A) 
was obtained. Since Ref. [17| assumed one-dimensional ge- 
ometry, the effective (Sharvin) resistance of the setup can 
be estimated as I/Gq. This demonstrates consistency of 
our results. It has been already argued in 18 that the 
rate of Ref. [l7] is suppressed by a factor accounting for 
the lateral spread of the ABS wavefunction in the leads. 
We have thus shown that the suppression factor is conve- 
niently expressed in terms of the effective resistance i?eff- 

Further, let us estimate the exchange coupling J as the 
interaction energy of two electrons (or holes) localized in 
the volume V. Owing to the strong screening in the bulk of 
the metal, the interaction quenches if the distance between 
the electrons exceeds the atomic scale a. The probability 
for electrons to be that close is estimated as a^/V. This 
yields J ~ i?at(a^/V), where we estimated the interaction 



by the atomic energy scale E^t- Since i' ~ (E'ata^)""'^ the 
estimation of J can be expressed in terms of the number 
N introduced above, 



J 



V 



A 

N 



A(i?effGQ) 



(9) 



We compare our estimations of F^^ and J to find that ex- 
change splitting exceeds by far the life-time broadening of 
the singlet state, hT^/J ~ A'^fiKei-ph — 10^**. We note 
that the ratio is independent of the effective resistance i?eff 
and is thus rather universal for all ABS setups. 

The last estimation we need is that of the SO split- 
ting ~ |eso|- In the case of a general superconducting 
constriction of size of the order of the scale of SO split- 
ting is given by a^^A [ij, where the SO coefficient in Al 
~ 10"^ can be interpreted as a probability to flip the 
electron spin in course of a scattering event. This however 
implies the energy dependence of the transmission ampli- 
tudes at a scale of ~ A. For a short break junction, this 
energy dependence is absent; this cancels the SO splitting 
^. The energy-dependent part of the transmission am- 
plitude comes about the scattering of an electron that has 
passed the constriction, proceeded at distance of the or- 
der of ^ in a lead and got back to the constriction (see 
Fig. [2]). The estimation for SO splitting is thus reduced 
by the backscattering probability. The latter can be again 
estimated as \/N ~ R^sGq. Thus, we find 



esol - «i?A/iV(A) ~ aiVAi^effGg 



so , 



(10) 



We compare spin and exchange splitting to find the for- 
mer to be a factor of lesol/-^ — — ^0"^ smaller not 
depending on the setup details. 

Let us provide typical scales of the estimated param- 
eters. For an experiment with i?off — 10 Ohm (that is, 
N ~ lO'^) and A//i ~ 10^^ Hz the above estimations yield 
J ~ 10® Hz, |eso| ^ 10^ Hz and F^ ~ 10"* Hz. We refer to 
these values in following concrete estimations. 

Spin-orbit relaxation of the triplet states. We 

proceed to describe the dynamics of the triplet TQTL 
states in conditions of time-independent magnetic field 
and superconducting phase drop at the junction. Quasi- 
particle annihilation in the triplet states is forbidden by 
spin blockade. However, the SO interaction given in eq. ([5| 
provides a relaxation mechanism by coupling the triplet 
states with the singlet. The coupling is described by the 
difference of SO pseudo- vectors €d- 

We compute the SO relaxation rate by treating the small 
quantities Cd/J — and hVd/J as perturbations. It is 
convenient to choose the spin quantization axis along the 
direction of vector {B+ts) and work with the parallel (per- 
pendicular) component of with respect to this direction, 

^d (^d )■ fi^'i decay rates Vt^ corresponding to the 
triplet states \Tj), j = — 1,0, -|-1. 

VT^,=Td{4f l{J±B)\ TT„=Td(ei) /J\ (11) 
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We estimate the decay rates — Td{a^)^ yielding 
— 1 Hz. We assume that the triplet states are long- 
lived at the scale of the measurement time. 

Singlet-to-triplet manipulation. We detail here 
the resonant manipulation of the singlet-to-triplet transi- 
tion. The pseudovector that defines the couplings be- 
tween the singlet and triplet states depends on the super- 
conducting phase if that can be modulated by an electric 
signal at frequency il, 



(f{t) = (fio + A [exp{iilt) + exp{~iiU)] 



(12) 



A being the dimensionless modulation amplitude A <^ 1. 
The corresponding modulation of the pseudovector is then 
given by: 

edit) = ediipo) + ^{'Pa)A[cxp{int)+cxp{~int)] . (13) 

The condition of resonant manipulation is achieved 
when hn matches the singlet-triplet energy spacing, that is 
different for different triplet states provided a sufficiently 
large magnetic field B |ego| is applied . In this case, 
each triplet state can be addressed separately. 

Let us find the change of probability to be in a triplet 
state as a result of a manipulation pulse of duration r. We 
employ rotating wave approximation (RWA) justified by 
B ^ |e'rf|,/ir(j. For the transition between |5) and \Tj); 
the resonant condition is = J + jB. It is convenient 

to introduce slowly varying states 
and 

equation 

i[Hs,p] - ^ [PST, 



Sj = exp(-iJt/2) IS") 
Tj) = exp[i(J/2 + jB)t] \Tj) to arrive at evolution 



S 



TdPss s) (S : hn, S 



h.c. 

+ h.c 



(14) 



where the manipulation amplitudes are = A(de^/ dip) , 
n± = A{de^/dip), e± EE -(±e^ + iey)/y/2. We assume 
that the manipulation starts at time t 3> T'^^ after the 
preparation, so that pss(O) = 0- If the initial probability 
to be in the triplet state is Pj (0) , the final one is given by 



(27^ - 1) cosh 



27/1 



1 sinh 



27 




-TdT - 1 



2(72 - 1) 



(15) 



The duration dependence is determined by the inverse di- 
mensionless strength of the pulse 7 = Td/iClj. For weak 
pulses 7 > 1, the singlet state decays faster than the 
coherent transition takes place and the triplet probabil- 
ity decays with the rate Fd(l — y^l — 7^-^) (Fig. [s^). At 
7 3> 1, this decay rate is twice the manipulation ampli- 
tude. For sufficiently strong pulses 7 < 1 one sees coher- 
ent oscillations of the probability (Fig.lsb) with frequency 



2f2j(l — 7^)^/^ on the background of the overall decay at 
rate F^. The full depletion of Pjir) can be achieved by 
tuning the pulse duration. 

Triplet-to-triplet manipulation. — Resonant ma- 
nipulation of the triplet-to-triplet transition is used to 
implement single qubit rotations. The coupling between 
|To) and cither of \T±) is realized both by SO interaction, 
6*3, and by the magnetic field, B, while the triplets \T^i) 
and do not mix. The manipulation can be achieved 

by modulating either coupling parameter. However, it is 
practical to modulate the magnetic field, as it can be much 
larger than the SO splitting in the junction. 

We consider a magnetic field with a static and an a.c. 
component oscillating at frequency 

B{t) = Bz + Bcxp{int) + B* exp(-im) (16) 

assuming B ^ B. The resonance condition is il — B/h, 
the same for both |T+i) and |T_i). In RWA, the effective 
Hamiltonian reads (B = B^ ~ iBy) 



Hs=\B\\To) {Ts\+h.c 
1 



IT, 



B \T+i) + B* \T^i) 



(17) 
(18) 



describing the Rabi oscillations with frequency Qji = 
\B\/fi between \To) and a superposition state \Ts). Af- 
ter a manipulation pulse of duration r the amplitudes ao 
and a, are transformed as 



ao(T) 




_as(T) 





cos(r2flt) — z sin(rij^i) 
-isin(riflt) cos(ilfl,i) 



ao(0) 



(19) 



Assuming the static magnetic splitting is of order B ~ J 
and the modulation is 10% of B, we estimate the Rabi 
frequency flu ~ 10^ Hz: this allows a fast and efficient 
triplct-to-triplet manipulation. The time available for the 
manipulation is set by the rate of the triplet decay Ft^ ~ 
IHz so that the number of rotations can be as large as 

Spin qubit. — Having understood the manipulation 
of the spin states, we can describe the design and operation 
of a simple spin qubit. The qubit states are |0) = |To) and 
|1) = \Tg). These together with the state defined below 
form an orthonormal basis of the triplet subspace. 



|2) 



V2\B\ 



B |r+i> - B* |r_i) 



(20) 



Let us give an example of an experiment with the qubit. 
In the beginning, the probabilities of all triplet states are 
the same, Po = Pi — P2 — T'/S, P being the probability to 
realize a long-lived triplet state introduced above. These 
probabilities can be manipulated by triplet-to-singlet tran- 
sitions. By tuning the pulse duration (Fig. |3|, we can 
achieve the full depletion of pa initializing the qubit to 
a mixed state with po = 0,pi = P2 = P/3. Subsequent 
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Fig. 3: Singlet-to-triplet manipulation. The upper plots 
show the pulse <^(t) of duration r; the insets show the time- 
dependence of the probability P to realize the long-lived triplet 
state. The lower plots depict the dependence P{j), in a) for 
the case A < Fd/ii^j, and in b) for the case A > Td/ii^j- 




Fig. 4: Qubit initialization, rotation and read-out. a) The 
three square pulses are applied in succession. The upper plot 
shows the two pulses of (p{t) used for initialization and read- 
out. The duration of the pulses is optimized to transfer all 
population from state |0) to \S}. The lower plot shows the 
magnetic field pulse of duration r used for triplet-to-triplet 
manipulation, b) The dependence P{t) of the probability to 
measure the lower critical current after the sequence of three 
pulses. 



rotation works if the qubit is in |1) bringing it to a super- 
position of states |0) and |1) by a triplet-to-triplet pulse of 
variable duration. Subsequently we read the resulting po 
after the rotation using the same pulse as used for the ini- 
tialization. Fig. |4] shows the result of applying the pulses 
for initialization, qubit rotation, and read-out in succes- 
sion. The measured probability to get the lower critical 
current P{t) = 2P/3 — Po{t), is plotted versus the pulse 
duration t manifesting Rabi oscillations. 

Conclusion. — We have presented a proposal of a 
novel qubit design using the spin states of two supercon- 
ducting quasiparticles trapped in a junction. Read-out of 
the qubit is based on spin-blockade that inhibits recom- 
bination of quasiparticles in the triplet state. We have 
described the resonant manipulation of singlet-to-triplet 
and triplet-to-triplet transitions and have explained the 
operation of the qubit. The qubit operation frequency 
estimated as £7 ~ 10^ Hz is much larger than the qubit 
relaxation rate F, il/F ~ 10^. Realization of our proposal 
would unambiguously demonstrate for the first time the 
spin properties of superconducting quasiparticles. 

* * * 

This work is part of the research program of the Sticht- 
ing FOM. 
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